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Solitons are non-dispersive wave solutions that 
arise in a diverse range of nonlinear systems, sta- 
blised by a focussing or defocussing nonlinear- 
ity. First observed in shallow water [Ij, solitons 
have subsequently been studied in many other 
fields including nonlinear optics, biophysics, as- 
trophysics, plasma and particle physics [2j. They 
are characterised by well localised wavepackets 
that maintain their initial shape and amplitude 
for all time, even following collisions with other 
solitons. Here we report the controlled forma- 
tion of bright solitary matter-waves, the 3D ana- 
log to solitons, from Bose-Einstein condensates of 
^^Rb and observe their propagation in an optical 
waveguide. These results pave the way for new 
experimental studies of bright solitary matter- 
wave dynamics to elucidate the wealth of existing 
theoretical work and to explore an array of poten- 
tial applications including novel interferometric 
devices [3j, the study of short-range atom-surface 
potentials |4] and the realisation of Schrodinger- 
cat states [5| ^6] . 

Bose-Einstein condensates formed from dilute atomic 
gases support bright soliton solutions for attractive inter- 
atomic interactions (focussing nonlinearity ) , manifesting 
themselves as localized humps in the field amplitude. In 
contrast, dark solitons appear as localized reductions in 
an otherwise uniform field amplitude, preserved by a de- 
focussing nonlinearity (repulsive interactions). The con- 
trol with which these systems can be manipulated, com- 
bined with the unique properties of matter- wave solitons, 
leads to a rich testing ground for theoretical descriptions 
of quantum many-body systems. Condensates are com- 
monly described by a mean-field treatment [7, 8 leading 
to the well-known Gross-Pit aevskii equation (GPE) in 
which the atomic interactions are described by a nonlin- 
ear term proportional to the 5-wave scattering length 
and the condensate density. In the one-dimensional (ID), 
homogeneous limit the GPE takes the form of a nonlin- 
ear Schrodinger equation which supports a spectrum of 
mathematically ideal soliton solutions. Experiments ap- 
proach this theoretically ideal scenario by confining the 
condensate in an elongated, prolate trap typically with 
tight radial confinement. However, this quasi-lD geom- 
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etry is usually accompanied by the presence of weak ax- 
ial harmonic trapping which removes the integrability of 
the system and prevents the appearance of true solitons. 
Nevertheless, solitary wave solutions remain which re- 
tain many similarities to the classical soliton solutions 
such as propagation without dispersion. 

Previously, bright solitary matter- waves have been re- 
alised in three separate experiments |1QHT2] . In each case 
a Feshbach resonance was used to switch the interactions 
from repulsive (a^ > 0) to attractive (a^ < 0) in order to 
form solitary waves out of the collapse instability [13 . In 
two of these experiments, multiple wavepackets were cre- 
ated, allowing the study of the dynamics during collisions 
in the trap. The observation of solitary waves raises many 
interesting questions regarding the relationship between 
the mathematical ideal and the experimental reality. It 
is unclear how soliton-like the solitary waves created in 
experiments with finite radial trapping and harmonic ax- 
ial confinement are. An answer to this question needs to 
be established before potential applications utilising soli- 
tary waves can be realised. At a more fundamental level 
it remains to be tested whether or not the GPE treat- 
ment fully describes the solitary waves created in exper- 
iments. Solitary waves realised experimentally typically 
contain ^1000 atoms, placing them well outside of the 
thermodynamic limit and potentially outside the reach 
of the mean-field description. Several theoretical studies 
of bright solitary waves beyond the mean-field descrip- 
tion have now been performed, either including effects of 
quantum noise using the truncated Wigner method [M] 
or using approximate analytic and numerical methods to 
simulate the full many-body problem [5l [15] . These gen- 
erate results potentially in confiict with the behaviour 
predicted by the GPE treatment. 

In this work we report the controlled formation of 
bright solitary matter- waves from a ^^Rb Bose-Einstein 
condensate. The experimental geometry is such that the 
velocity of the wavepackets can be precisely controlled, a 
key factor in facilitating the future exploration of solitary 
wave interactions and collisions. In addition, we observe 
and model the controlled refiection of solitary waves from 
a broad Gaussian potential barrier, demonstrating their 
particle-like nature. 

^^Rb is a prime candidate for solitary wave experi- 
ments owing to the existence of a broad Feshbach res- 
onance at ~155 G in collisions between atoms in the 
F = 2,mF = —2 state. We use this resonance to form 
a stable, repulsively interacting condensate in a crossed 
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FIG. 1: Expansion in the waveguide: (a) Experimental setup showing the glass science cell, the crossed dipole trap used to 
create the BEC, the optical waveguide and the quadrupole, bias and offset coils. Also shown in the cell is a super-polished Dove 
prism (blue), mounted on a macor support, to be used for future experiments, (b) Schematic of the release of the condensate 
from the crossed dipole trap into the waveguide, (c) Condensate expansion in the waveguide for as = 165 ao (black), 23 ao 
(red), 4 ao (blue), -7 ao (green) and -11 ao (purple). Solid lines are linear fits to the experimental data where the widths 
are rms values, (d) Condensate expansion rate in the waveguide as a function of atom number and scattering length. The 
solid line is the theoretical expansion rate calculated from a zero-temperature simulation of the experimental expansion using 
a cylindrically-symmetric, 3D Gross-Pitaevskii equation. As in the data, the expansion rate is defined using the change in the 
width of the BEC between 10 ms and 100 ms after release into the waveguide potential, which is approximately linear over this 
time interval in all cases 



optical dipole trap, shown in Fig.[TJa). The condensate is 
then loaded into a quasi- ID waveguide, better suited ge- 
ometrically to the observation of solitary waves. At the 
point of release into the waveguide, the magnetic bias 
field controlling the atomic scattering length is jumped 
to a new value (see Fig. [ijb)). As the BEC propagates 
along the waveguide the value of determines the rate 
of expansion of the condensate in the axial direction. We 
probe this expansion by measuring the condensate size as 
a function of time for different values of as shown in 
Fig. [ijc). Fitting the experimental data we can extract 
an expansion rate for the BEC, dependent on and N. 
This is shown in Fig.Illd), along with a 3D GPE simula- 
tion of the expansion (the solid line). At = — 11 ao and 
A/" = 2, 000 we see the expansion rate of the BEC becomes 
consistent with zero. This lack of dispersion with time 
indicates the formation of a bright solitary matter-wave. 

Figure [2] shows the propagation of this solitary wave, 
contrasted to that of a repulsively interacting BEC. As 



the repulsive wavepacket propagates the axial expansion 
causes a significant drop in optical depth not seen for the 
solitary wave. We observe a solitary wave propagating 
over a distance of 1.1 mm in a time of ~150 ms with 
very little distortion. 

To probe the stability of the solitary wave we investi- 
gate reflection of the wavepacket from a repulsive Gaus- 
sian barrier with a 1/e^ radius of 130 /im, shown in 
Fig. [3]^ a). Figures [sj^b) and (c) show the position of the 
solitary wave as a function of time in the presence of a 
760 nK barrier potential. In this case the barrier height 
is greater than the kinetic energy of the solitary wave and 
the wavepacket is cleanly reflected. 

Using a barrier much wider than the solitary wave size 
the atomic center-of-mass coordinate behaves classically, 
with the solitary wave acting as a single particle rolling 
up a potential hill. By varying the height of the potential 
barrier it is possible to select whether the solitary wave 
is reflected or allowed to travel over the barrier. The po- 
sition of the solitary wave after 150 ms is shown in Fig. 
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FIG. 2: Propagation in the waveguide: (a) As a repulsive BEG propagates along the waveguide the atomic interactions cause 
the condensate to spread, leading to a drop in optical depth, (b) In contrast, the attractive interactions present in a bright 
solitary matter- wave hold the atomic wavepacket together as it propagates, maintaining its shape with time. Grosscuts shown 
are the horizontal optical depth profiles of the condensates after 140 ms propagation time along the waveguide. 



[3[d) as a function of barrier height. The solid line is a 
theoretical trajectory, calculated using classical mechan- 
ics with no free parameters, showing excellent agreement 
with the data. 

In Fig. [sj^e) we compare the effect of reflection from 
the barrier for a solitary wave and a repulsive BEG and 
contrast the change in width to the case of a repulsive 
BEG propagating along the waveguide in the absence of 
the barrier. The solid lines are the theoretical condensate 
widths calculated by solving the 3D (cylindrically sym- 
metric) GPE. As expected, the solitary wave is robust 
against collisions with a repulsive Gaussian barrier and 
following the reflection maintains its shape, continuing 
to propagate without dispersion. In the absence of the 
barrier the repulsive BEG expands steadily in time. (We 
attribute the disagreement between experiment and the- 
ory at longer times to a small thermal component making 
the measurement of the condensate width less accurate.) 
In the barrier reflection case, an oscillation in the con- 
densate width is induced as a result of the larger spatial 
extent of the repulsive BEG causing it to be strongly com- 
pressed as it is reflected from the barrier. Such contrast 
in the behaviour of the repulsive BEG and the solitary 
wave reflection lends weight to previous theoretical pre- 
diction regarding the superior characteristics of solitary 
waves for observing quantum reflection from surfaces [4] . 

There is currently much theoretical interest [T6]-[20] in 
the scattering of solitary waves from narrow potential 
barriers where, if the barrier width is on the order of 
the solitary wave width, quantum effects are observable. 
At high kinetic energy soliton splitting is energetically 
allowed at narrow repulsive barriers. The effect of quan- 
tum tunnelling means the barrier can act as a beam split- 
ter p^, dividing the soliton into two parts [21^. These 



multiple wavepackets can then be used to investigate the 
phase dependence of binary collisions [17] , the behaviour 
of collisions of two solitary waves on a barrier [16l [18] 
and would provide a solid first step towards the reali- 
sation of a bright solitary wave interferometer. In the 
limit of low kinetic energy a mean- field GPE treatment of 
the problem begins to break down [22 and quantum be- 
haviour, described by the Lieb-Liniger Hamiltonian [23] . 
becomes more significant. Here, splitting of the soliton is 
energetically forbidden and it becomes possible to create 
Schrodinger-cat states [5\ [6] . 

The use of a narrow potential to controllably split a 
solitary wave presents an opportunity to investigate one 
of the key open questions arising from previous work; 
what governs the dynamics and stability of multiple soli- 
tary waves existing in the same trap? The long lived 
nature of the solitary waves and their apparent stability 
during binary collisions has been the subject of a wealth 
of theoretical work [14, 24 -27 . Within the framework of 
the GPE, the observed stability of soliton collisions can 
only be explained by imposing a relative phase (j) = n 
between neighbouring solitary waves [26] such that the 
collisions are effectively repulsive in character. Several 
other studies address the apparent stability of solitary 
waves in binary collisions, offering different interpreta- 
tions which do not require the imposition of a relative 
phase (j) = 7T between neighbouring solitary waves. The 
inclusion of quantum noise [14 or accounting for many 
body effects [TF both result in effectively repulsive in- 
teractions between solitary waves, irrespective of initial 
phase. Interestingly, incoherent, fragmented objects are 
also predicted to form in the many body formalism. Fur- 
ther experimental studies are undoubtedly required to 
address the role of the relative phase in solitary wave col- 
lisions and to test the different theoretical descriptions of 
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FIG. 3: Reflection from a repulsive Gaussian barrier: (a) Potential in the axial direction along the waveguide in the presence 
of the repulsive barrier. (Inset, upper: Combined waveguide and Gaussian barrier potential. Lower: Experimental setup.) (b) 
False colour images of a solitary wave reflecting from the barrier. The white line shows the location of the barrier centre, (c) 
Horizontal position of a solitary wave propagating in the waveguide in the absence (red) and presence (black) of the repulsive 
barrier, (d) The position of a solitary wave after 150 ms propagation time as a function of the barrier height. Red (black) 
points correspond to the solitary wave travelling over (being reflected from) the barrier. Solid lines in (c) and (d): Theoretical 
trajectory calculated using a classical particle model with no free parameters, (e) Condensate width following reflection from 
the barrier. In the absence of a barrier, a repulsive BEC will expand as it propagates (red). With the barrier in place, an 
oscillation in the condensate width is set up following the strong compression of the condensate at the barrier due to the shape 
of the potential (black) . A solitary wave undergoing the same collision emerges unaltered (blue) . Solid lines are the theoretical 
condensate widths calculated by solving the 3D (cylindrically symmetric) GPE. 



quantum many-body systems. 



METHODS 



Production of a tunable Bose-Einstein 
condensate 



Although reflection and splitting experiments show the 
potential to settle the theoretical debate over the solitary 
wave formation and dynamics, the ability to probe such 
narrow and hence rapidly varying potentials using these 
wavepackets also lends itself to an obvious application in 
precision measurement. Atoms close to a surface are sub- 
ject to the short-range Casimir-Polder and van der Waals 
potentials which can be measured using the classical and 
quantum reflection of bright solitary matter- waves [4 . 
Our apparatus includes a super-polished Dove prism for 
such studies, see Fig. [ija). Further in the future, the 
ability to deliver and manipulate ultracold atoms near to 
a solid surface may open up new routes to probe short 
range corrections to gravity [28 due to exotic forces be- 
yond the Standard model. 



We create a Bose-Einstein condensate with tunable 
atomic interactions using the method described in [29^. 
A magnetic Feshbach resonance is used to tune both the 
elastic and inelastic scattering properties of the atomic 
sample to achieve efficient evaporation. Importantly the 
resonance at 155 G in collisions between ^^Rb atoms in 
the F = 2,mF = — 2 state gives control over the s-wave 
scattering length close to the zero crossing of ~ 40 ao/G. 

The use of a magnetic Feshbach resonance means it 
is advantageous to work with a levitated crossed opti- 
cal dipole trap. This is formed from a single 10.1 W, 
A =1064 nm laser beam (IPG: YLR-15-1064-LP-SF) used 
in a bow-tie configuration as shown in Fig. [ija). The 
term 'levitated' refers to the use of an additional mag- 
netic quadrupole field whose vertical gradient is set to 
just less than that required to support atoms against 
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gravity. This trap allows the magnetic field, and hence 
scattering length, to be changed independently of the 
trapping frequencies. 



B. Loading the optical waveguide 

To investigate the creation of solitary waves we begin 
by forming a BEC containing up to 10,000 atoms at a 
scattering length of ^ 300 ao- The crossed beam trap 
in which the BEC is created has a roughly spherically 
symmetric geometry at the point of condensation, with 
final trap frequencies of ujx^y^z = 27r x (31,27,25) Hz. 
This trap is ill-suited to the observation of bright soli- 
tary matter-waves and thus we transfer the condensate 
into a more quasi- ID waveguide created by an additional 
1064 nm laser beam, focused to a waist of 117 fim and 
intersecting the crossed trap at 45° to each beam. This 
enters the glass science cell through the back surface of 
an anti-reflection coated fused silica Dove prism (to be 
later used for the study of atom-surface interactions [4 ). 

To load the condensate into the waveguide the scat- 
tering length is ramped close to = in 50 ms thus 
reducing the condensate size and creating a BEC ap- 
proximately in the harmonic oscillator ground state of 
the crossed trap. The BEC is then held for 10 ms to al- 
low the magnetic field to stabilise before simultaneously 
switching the waveguide beam on, the crossed beams off 
and jumping the quadrupole gradient in the vertical di- 
rection from = 21.5 G cm~^ to 26 G cm~^. Although 
it is advantageous in terms of the evaporation to be un- 
der levitated during the condensation phase, we must 
increase the gradient once we wish to transfer the atoms. 
This ensures a truer levitation of the atoms in the waveg- 
uide trap, thus maximising the trap depth of the beam. 
In addition, the presence of the quadrupole gradient pro- 
vides much of the, a lbeit weak, axial trapping along the 
beam, coaxial = 1/2^ jiW^ /ttiBq ^ 27rxl Hz [30 . Here, /i 
is the magnetic moment of atoms with mass m and is 
the magnetic bias field. The waveguide beam itself con- 
tributes < 0.1 Hz to the axial trapping, hence the mag- 
netic confinement dominates in this direction. At a beam 
power of 0.17 W the waveguide and quadrupole potential 
produce a trap of oOx^y^z = 27r x (1,27,27) Hz. Here the 
radial trap frequency {oOy^z) approximately matches that 
of the crossed beam trap at the point of condensation. 



C. Propagation in the waveguide 

A small offset (2.6 mm) between the crossed dipole 
trap, i.e. the waveguide loading position, and the 
quadrupole centre means that once loaded into the 
waveguide, the BEC propagates freely towards the mag- 
netic field minimum along the direction of the waveg- 
uide, undergoing harmonic motion. As the BEC propa- 
gates its rate of expansion in the axial direction is deter- 
mined by the scattering length. Although strictly speak- 



ing the expansion is non-linear over the full range of times 
measured, a linear approximation is valid over the range 
10 ms< t <100 ms from which we can extract a 'rate'. 



D. Control of the solitary wave velocity 

The position of the magnetic field zero in the axial di- 
rection of the waveguide can be displaced by an amount 
determined by the magnetic field gradient in this direc- 
tion, B' and a moderate offset field, ^offset, according 
to Ax = 5offset/(^V2) [31^. In this way the amplitude, 
and hence velocity, of the solitary wave motion can be 
precisely controlled due to the dominance of the magnetic 
potential over the optical confinement of the waveguide 
along the axial direction. The maximum velocity is given 
hy V = AcJaxiai where A is the amplitude of the motion, 
set by the separation between the minimum of the mag- 
netic potential along the axis of the waveguide and the 
release point from the crossed dipole trap. Using this 
technique the solitary wave can reach velocities of tens 
of mm s~^ when travelling through the centre of the har- 
monic potential or, alternatively, be brought to a near 
standstill, achieving velocities < 0.5 mm s~^. 



E. Classical reflection from a Gaussian barrier 

To produce the repulsive potential barrier we use a 
532 nm Gaussian laser beam (derived from a Laser Quan- 
tum Finesse laser), focussed to a waist of 131 /im hori- 
zontally and 495 /im vertically, with a power of up to 
2 W. The barrier is aligned to cross the waveguide in 
the horizontal plane at an angle of ~45° and is offset by 
455 /im from where the BEC is released from the crossed 
dipole trap, see Fig. [sja). This angle is restricted by the 
available optical access close to the trap centre. 



F. Theoretical modelling 

The release of the BEC into the waveguide poten- 
tial, and its subsequent expansion, was modelled at zero- 
temperature by solving the Gross-Pitaevskii equation in 
3D using a cylindrically-symmetric Fourier pseudospec- 
tral method. In all cases the initial non-interacting 
ground state of a harmonic trap with axial (radial) fre- 
quency 30 (27) Hz (corresponding to the crossed dipole 
trap potential) was released instantaneously into another 
harmonic trap with axial (radial) frequency 1 (27) Hz and 
offset by 2.6 mm along the axial direction (corresponding 
to the waveguide potential). The scattering length was 
instantaneously changed to the appropriate value of 
at the time of release. 

In cases where the barrier was present this was mod- 
elled as a Gaussian 'light-sheet' potential centered on a 
plane perpendicular to the axial direction, offset from 
the initial harmonic trap by 2.145 mm, and with height 
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760 nK and width 131 \/2 fim. Compared to the ex- 
perimental barrier beam this model neglects the vertical 
width of the beam which is large compared to the radial 
extent of the BEC in the waveguide, and includes the 
geometric factor a/2 to account for the 45 degree angle 
of the beam. 

Expansion rates were calculated from the full- width at 
half maximum of the BEC axial density profile predicted 
by the GPE (obtained by integrating over the radial coor- 
dinate) after 10 ms and 100 ms of expansion. In all cases, 
the change in radius over this time interval was approxi- 
mately linear. For the simulations in Fig.jsje), the width 
was calculated by convolving the BEC axial density pro- 
file predicted by the GPE with a 10 /im width Gaussian 
(to account for finite imaging resolution), and fitting a 
Gaussian distribution to the resulting profile using non- 



linear least-squares. 
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